Nodal Spin Density Wave and band topology of the FeAs based materials 



o 
o 

(N 



Ying Ran, 1,2 Fa Wang. 1 ' 2 Hui Zhai. 1 ' 2 Ashvin Vishwanath, 1, 2 and Dung-Hai Lee 1,2 

1 Department of Physics, University of California at Berkeley, Berkeley, CA 94720, USA 
2 Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 

(Dated: March 15, 2009) 

The recently discovered FeAs-based materials exhibit a (ir, 0) Spin Density Wave (SDW) in the 
undoped state, which gives way to superconductivity upon doping. Here we show that due to an 
interesting topological feature of the band structure, the SDW state cannot acquire a full gap. This 
is demonstrated within the SDW mean-field theory of both a simplified two band model and a more 
realistic 5-band model. The positions of the nodes are different in the two models and can be used 
to detected the validity of each model. 
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I. INTRODUCTION 

Since the discovery of superconductivity in 
Lai-zF^FeAs at T c = 2QK-, there has been mounting 
excitement associated both with the rapidly increasing 
T c when La is substituted by other lanthanoids^, as 
well as the similarity with the cuprate superconductors. 
As with the cuprates, the FeAs materials are quasi 
2D square lattice based transition metal compounds, 
which are magnetically ordered at stoichiomctry. On 
doping, in both cases, the magnetism is replaced by 
superconductivity. 

However, there are several significant differences. Most 
importantly, the FeAs stoichiometric compounds are not 
insulating^. Also, the magnetic order here is along the 
(it, 0) direction, and has a small moment^. In contrast to 
the cuprates, where only the d x 2_ y 2 orbital of the five Cu 
d-orbitals is important, the multi-orbital nature of the 
FeAs materials have been emphasized in several recent 
calculations^^. In this paper we point out a topological 
aspect of the band structure, closely connected with the 
multi-orbital nature of the material, that has important 
ramifications for the phases in this system. In particu- 
lar we show that a symmetry enforced band degeneracy 
at high symmetry points in the Brillouin zone leads to 
a band structure with nontrivial topology. This can be 
quantified in terms of a 'vorticity' quantum number. An 
example of vorticity ±1 is the Dirac node. Here the vor- 
ticity takes on values ±2. While such degenerate points 
occur in the band structure of other multi-orbital sys- 
tems, here they actually occur close to the Fermi level. 
In this article wc discuss both a simplified 2-band model 
and a realistic 5-band model where this is explicitly re- 
alized. 

An important consequence of this nontrivial band 
topology is that it leads to an unusual Spin Density Wave 
(SDW) state that is necessarily gapless. Specifically, the 
SDW wavevector connects hole pockets with vorticity 
= ±2 with electron pockets with zero vorticity. This mis- 
match forces nodes in the SDW gap function even in the 
presence of perfect nesting. Away from nesting the nodes 
are offset from the fermi energy, resulting in fermi pock- 



ets. We term such magnetic order 'nodal-SDW. We em- 
phasize that the topological feature required for the nodal 
SDW state exists both in other simplified models in the 
literature^ as well as first-principles band structures^. 

In the following we first derive and study a simplified 
two-orbital tight-binding model motivated by the quan- 
tum chemistry. Other studies£&i£ have also focused on 
two orbital models. However it has been argued^ that 
one needs at least three orbitals to accurately reproduce 
the LDA band structure and Fermi surfaces, e.g. location 
of the hole pockets in the Brillouin Zone. We then di- 
rectly study the more general and realistic 5-band model 
based on LDA calculation^. We consider the stoichiomet- 
ric compound (zero doping), and study the mean-field 
SDW phase. This mean field analysis confirms the exis- 
tence of nodes in the SDW gap function in both models . 
We also established the topological stability of the nodes. 
The location of the nodes in momentum space, and the 
associated fermi surface topologies, however, are differ- 
ent in the two models. This can be used to detect the 
validity of 2-band or 5-band model. A recent numerical 
renormalization group study by us has found precisely 
this nodal structure in the two-band model. The general 
requirements for the existence of these nodes are also dis- 
cussed. 



II. NODAL SDW IN THE 2-BAND MODEL 
A. 2-band microscopic model and band structure 

The Fe atoms form a square lattice whose principal 
axes are denoted as x and y, and the crystal structure 
axes arc labeled as X and ^[Fig. HJb)]. Let us first as- 
sume that Fe 3d xz and 3d yz are the relevant orbitals to 
describe the low energy physics of this material. Because 
of the tetragonal symmetry they are locally degenerate, 
and we use their linear combinations 3dxz and Zdyz as 
our basis since they have clear symmetry when hybridized 
with the nearest As 4p orbitals. Other d orbitals will be 
ignored here, but included in a later section. 

In a simple chemistry picture[Fig. [IJa)] all symmetry 
allowed Fe 3d- As 4p hybridizations are assumed to domi- 



2 




FIG. 1: (color online) (a) The 3dxz and 3dyz orbitals of the 
Fe atoms and the 4pz orbital of the As atom. The blue(red) 
arrows represents one second order process contributing to 
t2(ti). (b) The tight-binding hoppings between the d-orbitals. 

nate over direct Fe 3d hybridization. This naturally leads 
to large nearest- neighbor (NN) orbital-changing hopping 
fi (see Fig.QJb), note that t\ has opposite signs between 
vertical and horizontal bonds), next-nearest-neighbor 
orbital-preserving hoppings t% and t' 2 . If we further as- 
sume that the hybridizations between 3dxz,YZ an d 4p^, 
shown in Fig. [lja) dominate, we expect £2 ~ ti > l^li 
while the direct hopping t[ is expected to be much smaller 
than these three^. Given the empirically observed SDW 
order, we take — t' 2 > \ti\, which leads to nested elec- 
tron and hole pockets at (tt, 0). 

To simplify the notation we will use d\ and d,2 for elec- 
tron operators associated with the XZ and YZ orbitals. 
The tight binding Hamiltonian has only one Fe atom per 
unit cell, and we choose the Brillouin zone accordingly. 
Many other studies work with a two Fe atom unit cell, 
a comparison requires an appropriate folding of the Bril- 
louin zone of the present study. Its Fourier transform is 
(independent of spin) 

k 

where the sum is over k x 6 [— tt, tt), k y £ [— 7r,7r), the 
2x2 matrix K(k) is 

K(k x , k y ) = 2t\(cosk x ~ cosk y )T\ — 2(t2 — t 2 ) sinA^ sin 
+ [2(f 2 + t' 2 ) cos k x cos k y + 2t\ (cos k x + cos k y )\ • 1 (2) 

and 71,2,3 are Pauli matrices. As argued previously we 
expect t<i > t\ 3> t 2 ,t[ > 0. The two energy eigenvalues 
of Eq.© are: 

E± (k) = 2(t% + t 2 ) cos k x cos k y + 2t\ (cos k x + cos k y ) 
± 2^j t\ [cos k x - cos ky] 2 + (t 2 - t' 2 ) 2 sin 2 k x sin 2 k y (3) 

In Fig©a) we plot the band structure Eq.Q. At half- 
filling the fcrmi level will cut out two hole-like fermi sur- 
faces around (0, 0) and (tt, tt), and two electron- like fermi 
surfaces around (tt, 0) and (0, tt) (see Figl^b)). Note the 
band touchings at (0, 0) and (tt, tt). This endows the hole 
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FIG. 2: (color online) For ti = l,t 2 = 1.7,4 = 0.3, ti = we 
plot (a): the dispersion E±(k) as shown in Eq.Q together 
with the half-filling fermi energy Ef , and (b): the perfectly 
nested fermi surfaces as given by Eq.@ and the direction of 
vector (cos <f>%, sin </>g) defined by rewriting K(k x , k y ) in Eq.© 
as K(k x , k y ) — a(k)l + b(k) [sin </>jTi +cos 0jts] where b(k) > 
0. 

Fermi surfaces with 'vorticity' ±2, where the spinor de- 
scribing the admixture of 3dxz and 3dyz orbitals rotates 
twice on encircling these Fermi surfaces. The simultane- 
ous presence of inversion and time reversal symmetry in 
this band structure allows us to choose at each k point, 
real spinor wavefunctions which are hence confined to 
a plane. Vorticity in this spinor field is therefore topo- 
logically protected - the singularity at the vortex center 
forces the orbital degeneracy at (0, 0) and (tt, tt). In con- 
trast, the electron Fermi surfaces are topologically trivial, 
with no winding as shown in Figure^. This topological 
characterization of the fermi pockets is also present in 
more realistic LDA calculations^^, although it has not 
been previously commented upon. 

At half filling, the total electron pocket area equals 
the total hole pocket area. If the pockets are every 
small, i.e. t 2 w t\, Taylor expansion of Eq.Q gives four 
nearly circular fermi pockets with same area. Therefore 
the electron and hole pockets are nested at momentum 

MVCM). 

It is interesting to consider for a moment the case 
t[ =0. Then, the model decouples into two independent 
.t\ — t 2 — t 2 checkerboard models. Each give rise to one 
liole and one electron pocket, separated by (tt, tt) from 
the other electron-hole pocket pair. Unexpectedly, the 
electron and hole pockets are precisely nested with mo- 
mentum (tt, 0)/(0, tt), as long as t\ = 0, even when the 
pockets are large and non-circular. One can check that 

2t^ 

at the half filling Fermi energy Ep = u _^ t , , the fermi sur- 
face wavevectors k F of the two bands: £L (k F ) = Ep and 
E+(k F + (tt, 0)) = Ef satisfy exactly the same condition: 

[(f 2 + t' 2 f cos 2 fcf - t\] [(t 2 + t' 2 f cos 2 k F y - t\\ 
=(*a + t' 2 f(t 2 - t' 2 ) 2 sin 2 k F sin 2 k F . (4) 

When fj > the two hole pockets have different sizes. 
The electron pocket around (tt, 0) ((0, tt)) is elongated 
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along the k y (k x ) direction and the perfect nesting is 
lost. 



B. Mean Field Study of the SDW order in the 
2-band model 



Now we include on-site interactions in an extended 
two-band Hubbard model with H = Hq + Hi where, 
as in Ref@: 

Hi = - Y^( n ii + n h) + (U- 2 J) n «i^2 

i i 
i 

+ J 5Z ( d Ii,T d ii,l d i2a d i2,T + h.c.) (5) 



where the first and second terms are the intra-orbital 
and inter-orbital Coulomb repulsions. The third term the 
Hund's coupling and the fourth term is the inter-orbital 
pair hopping. First principles calculations^ suggest that 
the FeAs material is in the intermediate coupling regime 
and we choose t = 1, t[ = 0.2, i 2 = 1.7, t' 2 = 0.3 for 
our Hopping Hamiltonian Hq, where we estimate t\ ~ 
0.3eU to get the right bandwidth. Also, we use U = 4 
(1.2eV), J — 0.4, a factor of three smaller than in£ to get 
reasonable SDW transition temperatures and moments 
compatible with experiments. 

For small t' 1; where good nesting prevails and in the 
presence of repulsive onsite interactions, it is natural 
to consider SDW order at wave- vector (tt,0). Then, 

H S dw = MabYsii-Y* ( d lat d iM ~ d [ a i d i,bl) where the 
spin direction is assumed to be along the S z axis. Be- 
cause of the multi orbital nature of the system, differ- 
ent flavors of SDW are allowed, described by the Hermi- 
tian matrix M a b, which may be parameterized by four 
real numbers M ab = [4> t + 0iTi + 2 t 2 + <for 3 ] a b where 
(To)ab = 5 a b- We perform the finite temperature mean- 
field study by using a trial density matrix of the mean- 
field Hamiltonian Hmf = Hq + Hsdw- We construct a 
trial free energy based on this mean-field density matrix: 
ftriai(0o, 01,02,0a) = F M f + (Hi - H S dw)mf where 
Fmf is the free energy of the free fermion system de- 
scribed by Hmf- The Feynman inequality 12 F < F tr iah 
implies that we need to minimize F tr i a i over the mean 
field parameters <pi keeping the electron density fixed. 

Implementing this we find that the model Hq + Hi has 
a unique SDW phase at low temperature characterized 
by 0o 7^ 0, 0i 7^ 0. The symmetry of this phase is con- 
sistent with the regular (w, 0) magnetic order depicted 
in Fig|4] If we denote the order parameter operators 
m ab 



= d 



d iai d h b l> m FigI3 a ) we plot these 



magnetizations as a function of temperature and it is 
clear that (mn) = (7^22} 3> (^12)- The fact that the 4>2 
and 03 orders are not mixed in can be understood from 
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FIG. 3: (color online) For ti = 1, A = 0.2, i 2 = 1.7, t' 2 = 
0.3, U = 4 and J = 0.4 and half-filling, we show (a): the 
magnetization (mil) = (^22) and (m.12) of the SDW phase 
as a function of temperature. The mean-field calculation is 
done on 20 by 20 lattice. (b): the zero temperature SDW band 
structure. We plot the dispersion of two lowest energy bands 
in the reduced Brillouin zone, (c): the fermi pockets. The 
four nodes are not at the same energy and therefore the half- 
filled system has small fermi pockets. The two nodes close to 
k y = are hole doped (blue pockets) and the two nodes close 
to k y = 7r are electron doped (red pockets) . 



symmetry Under reflections P y about the y-axis cross- 
ing the plaquctte center, 0o and 4>i transform differently 
from 03. Under time reversal, 02 order transforms dif- 
ferently from a SDW, and in fact describes a spin-orbital 
locked state. 

The band structure resulting from this SDW is shown 
in FigJ3£d), where electron and hole pockets in the same 
direction as the SDW axis are visible for each pair of 
Fermi surfaces. Even with perfect nesting, the Fermi 
surfaces arc not fully gapped, instead there are Dirac 
nodes as in FigJ3Jb). In the absence of nesting these 
Dirac nodes are at slightly different energies and the half- 
filled system has small fermi pockets. As shown below, 
the presence of such pockets are required under fairly 
general conditions . This is the main point of the current 
study. The presence of such pockets are detectable by 
ARPES, and also, due to the velocity anisotropy present 
for the nodes, via conductivity measurements. Based on 
Drudc's formula we compute the DC conductivities in x 
and y directions for a mean-field SDW state with nodes 
along the x direction and with a T = magnetic moment 
of 0.3/xb per Fe atom4. We find a xx /a yy = 6.2. Any state 
with broken rotation symmetry would have conductivity 
anisotropy, but the large value here stems from the nodal 
structure discussed. 

For the other SDW type orders, 02 or 03, the nodes 
are along the direction orthogonal to the SDW axis. Note 
that although the nodal structure of our SDW gap func- 
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tion resembles that of a p-wave symmetry, the SDW order 
that we found is completely on-site and inversion sym- 
metric. 



C. No-full-gap 'theorem' in the two band model 

In the following we argue that the Dirac nodes are 
topologically stable as long as the SDW order satisfies 
three conditions (a brief account of this argument was 
presented in Refpii): (1) collinear order(denote the mag- 
netization direction by h) (2) inversion (about the Fe 
site I) symmetry and (3) Effective time reversal symme- 
try T1Z = S7Z(h — > — n) o TTZ obtained by combining 
time reversal and spin reversal (T1Z is time-reversal and 
SlZ(h — > —n) is the 180° spin rotation which flips the 
direction of magnetization). These three conditions are 
naturally satisfied by a (tt, 0) collinear SDW which is con- 
sistent with experiments*, and the mean-field SDW that 
we find also satisfies these conditions. 

Since we focus on the consequences of the nontriv- 
ial band structure, we begin by turning on a very weak 

SDW order: M = T,k a U Mab ^ d l a ,k d b/3,U^,o) ( where 
we have rotated the magnetization to the S z direction). 
We therefore need to consider degenerate perturbation 
after folding the Brillouin zone by (tt, 0). The band 
crossing between the two bands defined in equation [3j 
i.e. E-(k) and E + (k + (tt, 0)) forms a loop in momen- 
tum space. In the perfect nesting case these loops are 
identical to the electron or hole fermi surfaces. Let us 
focus on the loop around (0,0), and denote the mo- 
menta on this loop by k*: thus E_(k*) = E + (k* + 
(tt, 0)). Let us call the hole pocket wave-functions on 
this loop \iph(k*) = \ip-(h*)) and the electron pocket 
wave-function \ip e (k*) = \tp+(k* + (tt, 0))). Now the de- 
generacy on this loop is lifted when the matrix elements 
of the SDW order between these two kinds of states: 
m(k*) = (iph{k*)\M\il) e (k*)) is non zero. 

Now, I symmetry requires M a t(—k) = M a b(k). Also, 
given the winding of the hole Fermi surface wavefunction 
shown in Figure^, under inversion we have \if>h(~ k*) = 
-\tp h (k*)} while \ip e (-k*) = \tp e (k*)}. Putting this to- 
gether we have m(— k*) = —m(k). In addition, T1Z 
symmetry requires M a b(k) to be real hence so is m(k*). 
We thus conclude that m(k*) must have at least two 
sign changing points, K and —K, on the band crossing 
loop. These are the two of the Dirac nodes in the SDW. 
Similarly there are another two Dirac nodes on the band 
crossing loop around (0,7r). 

Generically the nodes are not at the chemical potential, 
and leads to Fermi pockets. However, a Fermi pocket 
deriving from a node is known, e.g. from the context of 
graphene, to be different from a regular Fermi pocket. In 
particular, electrons acquire a nontrivial Berry's phase of 
tt on circling such a nodal pocket. We have thus argued 
for the stability of the Dirac nodes in SDW based on 



the two-band model Eq. ([2]) . We now show that a similar 
result holds for the more realistic five band model, again 
the topology of the quadratic bands that touch at the T 
point are responsible for this result. 



III. NODAL SDW IN THE 5-BAND MODEL 

While the two band model serves as a useful guide 
to the nontrivial physics in the SDW state, it differs 
from LDA calculations of the electronic structure of these 
materials^ in important ways. This is easiest to see in the 
unfolded Brillouin zone scheme, with a single Fe atom per 
unit cell (here unit translations along the x and y axes are 
followed by reflections in the xy planed). While the two 
band model has a hole pocket at T — (0, 0) and another 
one at (tt, tt), the LDA calculation predicts two hole pock- 
ets around the T point. The electron pockets although 
centered around the same locations in both cases, acquire 
a d xz + d xy character at (0, tt) and a d yz + d xy character 
at (tt,0) in the LDA calculations. Hence they are also 
rather different from the two band model that does not 
include the d xy orbitals. A 5-band hopping Hamiltonian 
including all the iron d-orbitals is required to capture the 
fermi surface topology of the LDA calculation. 

In this section we will show that even in the five band 
model, the SDW is necessarily gapless (nodal-SDW), de- 
spite these important differences. An important role here 
is played by the fact that the two hole pockets at the T 
point are derived from the d xz , d yz , which are precisely 
the orbitals that enter the two band model. In this case 
we will prove that there must be at least two Dirac nodes 
close to fermi level in the SDW phase. 



A. The mean-field study of SDW in the 5-band 
model 

We again apply the trial density matrix method to 
study the (tt, 0) SDW instabilities of the 5-band model. 
We take the 5-band hopping Hamiltonian Hq from 
Kurokiei al£- (Eq.(l) and Table I). We then turn on a 
on-site interaction in the following forir& 

H i = U 2J niafUiai + (U - 2 J) ^2 n la n ib 

i.a i,a<b 
a,/3=T4 

i,a<ib 

+ J Y1 (4a,i d L,i d *bA d ibA + h.c) (6) 

i,a<b 

where a, b = 1,2... 5 label the 5 orbitals 
{d 3Z 2-R2,dxz,dYZ,d X 2~Y2,dxY}- We have pa- 
rameterized the inter-orbital Coulomb interaction by 
U — J. While this is strictly expected to hold within the 
t2 g and e g levels^, to reduce the number of interaction 
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FIG. 4: The real space pattern of the (tt, 0) SDW on the 
Fe-atom square lattice. 



parameters, we assume it for the five band model as 
well. 

We only consider the on-site SDW order, because the 
interactions taken to be on-site. Then, the (tt, 0) SDW 
order parameter with spins assumed to point along the z 
axis induces the following mean field term in the Hamil- 
tonian: 



H.SDW 



5 

E 

a,b=l 



M ab (4,aa^a'd lM ') (7) 



which is parameterized by M a b a Hermitian matrix with 
25 real parameters. The orbital structure can leads to 
many different SDW states which break symmetry in dif- 
ferent ways. Energetically, we find that the preferred 
state always has the same symmetry as the regular (tt, 0) 
SDW shown in FigQ] That is, the state breaks time re- 
versal (TTZ) and 180° spin rotation about an axis perpen- 
dicular to the ordering direction S1Z, but preserves their 
combination (TTZ'). Also the state is invariant under in- 
version (X) which is a 180° rotation in the x — y plane, 
the P x reflection around the x-axis crossing an Fe atom 
(which actually combines with z — > — z reflection and is 
a 3-D 180° rotation around the x-axis), and similarly the 
P y reflection (P x o P y = X) . At the end of this subsection 
we comment on other possible states, that break different 
symmetries. 

We now discuss details of the mean field solution, as 
the on-site interaction strength is varied. Since this is not 
accurately known, we note the resulting ordered moment 
and Fermi surface topology in each case, which may be 
directly compared with experiments. The Hund's cou- 
pling J is assumed to be about 20% of U. 

(1) U=1.0eV, J=0.2cV. The temperature evolution 
of the net SDW magnetization M = (X^o ^lo T^w.T — 
d\ a i dj ,l) (assuming g = 2 this is in unit of Bohr mag- 
neton) from mean field theory is shown in FigJSJa). 
The mean field transition temperature obtained is T c = 
0.026eV while the zero temperature magnetic moment 
obtained is ~ 0.23yUs. The latter is consistent with neu- 
tron scattering experiments on LaOFeAsX. In Figj5jc) 
we present the fermi surfaces of the zero temperature 
half-filled SDW phase with these parameters. Note that 
the double degenerate point at (0,0) of the 5-band hop- 
ping Hamiltonian is split into two Dirac nodes, C and 
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FIG. 5: (color online) (a) Magnetization with U = leV and 
J = 0.2eV, as a function of temperature (computation is per- 
formed on a 40 by 40 lattice with periodic boundary.) (b) 
T=0 dispersion in this SDW state showing the two low en- 
ergy bands along k x and k y axes around (0, 0). There are four 
nodes close to fermi energy: A and A' on the k x axis and B 
and B' on the k y axis. The double degenerate point at (0, 0) of 
the 5-band hopping Hamiltonian is split into two Dirac nodes, 
C and C , on the k y axis. We also plot the zero temperature 
fermi surfaces (blue pockets are hole-like and red pockets are 
electron-like) with (c) U = leV and J = 0.2eV and with (d) 
U = 1.2eV and J = 0.25eV and with (e) U = 1.4eV and 
J = 0.3eV. 



C", on the k y axis. The electron pocket at (0, tt) still 
intersect with fermi level and contributes significantly to 
the density of states. However, on increasing the interac- 
tion strength this feature is suppressed. Stronger repul- 
sive interactions will tend to gap out more of the fermi 
surface and the electron pockets at (0, tt) can be fully 
gapped out. We note that the area occupied by electron 
pocket is 3.5% of the magnetic Brillouin zone (defined by 
— 7r/2 < k x < ir/2 and — it < k y < tt). The hole pockets, 
of course, occupy the same area. 

(2) U = 1.2eV and J = 0.25eV. In this case we obtain 
a somewhat larger low temperature moment ~ 1.04/is. 
In FigO^d) we plot the fermi surfaces of the zero temper- 
ature SDW. As compared to the previous cases, we see 
that the fermi surfaces around (0, tt) disappears and the 
fermi level adjusts itself to form one hole pocket around 
(0, 0) and four electron pockets. The hole pocket occu- 
pies 2.2% of the magnetic Brillouin zone while the area 
of the electron pockets on the k x axis is 0.1% each, and 
of those on the k y axis is 1.0% each. 

The electron pockets on k x axis and on k y axis are 
fundamentally different in that the A^-pockets arise from 
the Dirac nodes A and A' below them, and thus pro- 
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tected. On the other hand the /c^-pockets can be easily 
gapped out by turning on interactions because they are 
simple band bottoms (after B(B') and C(C") annihilate 
each other, which already happens in the current case). 

(4) U = 1.4eV and J = 0.3eV. The zero temperature 
moment is now large ~ 2.3/jb. Now, the two electron 
pockets along the k y axis are completely gapped out and 
we only have a small hole pocket around (0, 0), and two 
electron pockets on the k x axis nearby (see FigJSfe)). The 
fermi surface topology and large moment obtained in this 
case are the closest to the LDA results^. The hole pockets 
occupy 1.4% of the magnetic Brillouin zone, while the 
electron pockets occupy 0.7% each. 

Note, for the larger interaction strengths, the area oc- 
cupied by the residual Fermi surface in the SDW state 
is very small, typically a few percent. In general, in- 
teractions that drive the SDW formation would tend to 
lower this area. One may naturally expect that this 
would rapidly lead to a fully gapped state on increas- 
ing U. However, as explained in detail in the next sec- 
tion, there is an intrinsic mechanism that blocks such 
a fully gapped state. A combination of symmetry and 
band topology necessarily leads to a gapless SDW state, 
over a wide range of coupling strengths. This provides a 
'natural' protection of the small pockets that appear here 
and in LDA calculations, which have now been observed 
in magnetic oscillation experiments^. 

Other Possible Orders: Before we conclude the mean 
field study of SDW orders, we comment on other kinds of 
(0, 7r) orders that could be stabilized with onsite interac- 
tions. The 25 parameter M a b matrix admits a plethora 
of different orders, which may be separated into the fol- 
lowing four classes, according to the symmetries of the 
resulting SDW Hamiltonian Hsdw'- 

(i) TTZ' even, X even and P x ,P y even. (6 parameters) 
This is the SDW state that was considered above. 
The no-full-gap theorem discussed below, applies 
to this case. 

(ii) TTZ' even, X even and P x , P y odd. (3 parameters) 

(iii) TTZ' even and X odd. (6 parameters) 

(iv) TTZ' odd. (10 parameters) This case is rather exotic 
because T1Z is even but STZ is odd. This is similar 
to the symmetry of a spin-hall insulatoi— . 

To compare the relative stabilities of these different 
states, we choose U = leV and J = 0.2eV and first per- 
form an unbiased minimization of all the 25 parameters. 
We find that class (i) is always the low free energy so- 
lution and with the highest T c = 0.026cV. Even if wc 
suppress order parameter (i) by hand, we find the sys- 
tem has no instability towards (ii), (iii) and (iv) down to 
O.OOOleV. (Computations were performed on a 40 by 40 
lattice with periodic boundary.) We conclude that SDW 
(i) is the low free energy phase in the model Eq.© and 
is consistent with the ordered pattern observed in exper- 
iments. Hence, we do not pursue studying these other 
kinds of SDW order. 



B. No-full-gap 'theorem' in the 5-band model 

In this section we explain why the nodal SDW found 
in the mean field study appears. We assume that the 
SDW has the symmetries in (i) above, as found in mean 
field theory. Briefly, we use reflection symmetry along 
the k x = and k y = to label bands with a reflection 
eigenvalue. Bands connected to the electron pocket and 
a hole pocket are forced to have opposite eigenvalues and 
hence do not split in the SDW state along these lines, 
leading to a gapless state. An important role is played 
by the band touching at the k = (0, 0) point, as in the 
two band model. While this reasoning holds for weak 
SDWs, we extend it to include strong SDW instabilities, 
where the location of the band intersections can migrate 
to the r point. Even in this case so we show that at least 
two gapless Dirac nodes will remain. 

Weak SDW Limit: The reflections P x and P y around 
the x and y axes passing through the Fe atoms 
(which actually combine with z — > — z reflection and 
are 3-D 180° rotations), act on the orbital basis of 
{d 3 z 2 -R 2 i dxz, dyz, d X 2-Y 2 ? dx y} as follows (note that 
our definition of X and Y axes given in Fig [T] is different 
from the definition in Kuroki, et aZ.— by a 90° rotation.), 
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Spin is left invariant as we ignore spin-orbit interactions. 
We first consider the P x reflection symmetry. Along k x 
axis the Bloch Hamiltonian is P x symmetric and the 
wave-function should be eigenstates of P x with eigenval- 
ues ±1. We thus can simply present the eigenvalue of 
P x of each band as shown in Figl6l If we focus on the 
three bands close to fermi level, we find that the electron 
pocket is P x odd, the large hole pocket is P x even and 
the small hole pocket is P x odd. Since the relevant SDW 
orders are P x even, the SDW induced gap along k x axis 
between the electron pocket and the large hole pocket 
must vanish, i.e. there is a band crossing. These are la- 
beled A and A' in the example of Fig[5jD. We note that 
the bands corresponding to the large hole pocket and the 
small hole pocket must have opposite P x eigenvalues and 
as result the nodes must exist no matter whether the 
electron pocket is P x even or odd. The simplest way to 
understand this is to note the double degenerate wave- 
functions at (0,0) are nothing but d xz and d yz . Hence, 
these band touchings play a crucial role here, as in the 
two band model. 

Similarly there also must be nodes along the k y direc- 
tion by studying the P y eigenvalues along the k y axis. We 
find that the electron pocket is P y odd, while the large 
hole pocket is P y even and small hole pocket is P y odd. 
Therefore we expect the SDW gap vanishes along k y di- 
rection between the electron pocket and the large hole 
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FIG. 6: (color online) We present the P x eigenvalues (blue for 
+1 and red for —1) of the 5 bands along k x axis. 



pocket as well. These are labeled B and B 1 in the exam- 
ple of FigO). Hence, in the limit of a weak SDW, the 
Fermi energy will inevitably cross at least one of the pro- 
tected bands, and a full gap cannot result. Note, the ex- 
istence of pockets is independent of how good the nesting 
is, in contrast with conventional SDWs which are gapless 
only due to the absence of nesting. A further distinction 
is that some of the gapless pockets in our case result from 
Dirac nodes. 

Strong SDW: The above four band crossings are ob- 
tained assuming the SDW is a weak perturbation. If 
SDW is strong can these nodes annihilate and lead to 
a fully gapped SDW? Note, in the presence of TTV and 
X symmetry, a band crossing can only be removed via 
annihilation with a partner. Since the band crossings 
appear close to the T point, it is particularly important 
to address whether they can annihilate by coming to- 
gether at that point. For completeness, it is important 
to note that in addition to the four band crossings dis- 
cussed above, there are two additional ones that arise 
when Cn v is broken in the SDW state. The double vor- 
tex at (0, 0) is split into two Dirac nodes with the same 
chirality when the SDW order is turned on. Therefore 
there are totally 6 Band crossings around (0,0) in the 
weak SDW limit (labeled as A,B,C and A',B',C in the 
mean field example of FigJSJx ) . In the Appendix we prove 
an important theorem based on a topological argument 
ensuring these 6 nodes can never annihilate in pairs at 
(0,0). We show that this is ensured by the fact that 
the two wavefunctions at (0,0), which corresponding to 
dxzidyz, are both odd under inversion (X). The nodes 
can only annihilate in sets of four. We conclude that 
there must be at least two Dirac nodes left. This result 
is topologically stable. Indeed, on going to stronger in- 
teraction strengths this is the state realized in our mean 
field study e.g. in FigO- In that case the two left-over 
Dirac nodes are along the k x axis. 

Now we compare the nodal SDW in 2-band model and 
5-band model. One main difference is that the fermi sur- 
face topologies and shapes are quite different. In par- 
ticular in the 5-band model there may not be a large 
anisotropy of conductivity as in the 2-band model be- 
cause at least there is a rather circular hole pocket around 
(0,0). The positions of Dirac nodes and fermi surface 



topologies, which may be easily measured by single crys- 
tal ARPES, can serve as a way to directly detect the 
validity of 2-band model or 5-band model. 



IV. CONCLUSIONS 

In this paper we studied the spin density wave (SDW) 
ground state of the undoped FeAs compound. We find 
that the combination of physical symmetry and the topol- 
ogy of the band structure, naturally stabilizes a gapless 
SDW ground state with Dirac nodes. We first study a 
popular two-band model due to its simplicity, where this 
mechanism is manifest. We also study the more realistic 
five-band model, where the same result obtains. These 
two rather different models share a key topological fea- 
ture of the band structure: the double degeneracy at 
k = (0, 0) enforces a wavefunction winding around this 
point in the Brillouin zone. 

In both models we perform the mean-field study al- 
lowing for all possible collincar magnetic orders at (w, 0), 
and find the inversion and reflection parities of the low- 
est energy magnetic ordered phase. We then show that 
the SDW ground state in both models have stable Dirac 
nodes protected by the inversion symmetry and the 
topology of the band structure. These Dirac nodes are 
close to Fermi level and thus may be directly observ- 
able in ARPES experiments and might also control the 
low energy thermodynamic and transport properties of 
compound. They arise due to the vanishing of the SDW 
matrix elements along a high symmetry line in the Bril- 
louin zone, which leaves the Fermi surfaces ungapped in 
this direction. We also proved a general result on the 
stability of Dirac nodes against pairwise annihilation in 
an inversion symmetric system (Appendix [A]) which may 
be applied to more general situations. While strong in- 
teractions tend to increase the SDW gap and reduce the 
Fermi pocket area, the nodal nature of the SDW does 
not allow a full gap to open over a wide range of interac- 
tion strengths. Hence, one expects to be left with small 
residual Fermi surface pockets, which naturally explains 
the small Fermi surface areas (0.52% and 1.38% of the 
Brillouin zone with two Fe atoms per unit cell) observed 
in magnetic oscillation experiments^. 

Although we find stable Dirac nodes in the SDW 
ground state in both the two and five band models, the 
number and the locations of the fermi pockets are dif- 
ferent. These differences can serve as ways to determine 
which is a better model of the material. Effective low 
energy theories of the FeAs materials should ideally in- 
corporate the nodal nature of the SDW state, which may 
also have important consequences for other phases in this 
system. 
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an early stage and support from LBNL DOE-504108, 
Insf-dmr 06456911 D. H. L. was supported by DOE DE- 
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APPENDIX A: ANNIHILATION CONDITION 
FOR DIRAC NODES 

We consider the general question of when a pair of 
Dirac nodes can come together and annihilate, to give 
rise to a non-singular band structure. We assume the 
existence of both time reversal and inversion symmetry. 
We consider the case where a pair of Dirac nodes come 
together at a point M in momentum space - which is 
constrained by the symmetries above to be invariant un- 
der inversion. When the Dirac nodes are near this point, 
we can restrict attention to just the two bands that make 
up these nodes. At the M. point, they can be labeled by 
their eigenvalues under inversion I\ = ±1 and I2 = ±1. 
Below we show that only if /1/2 = — 1 c & n the pair of 
nodes annihilate. Otherwise, they necessarily lead to a 
band touching with quadratic dispersion at the M. point, 
when they are brought together. 

To derive this result, we use a Berry's phase formula 
to fix the inversion eigenvalue of the Ai point states be- 
fore and after node annihilation. This places the required 
constrain on the inversion eigenvalues, if the nodes are to 
annihilate. We first define the inversion parity o~{C) of a 
half loop C connecting two points P and P' , which are 
mapped to one another by inversion. We start with the 
wavefunction of one of the bands at P, "Jp), which can 
be taken to be real given that we have both time reversal 
and inversion symmetry. This is evolved adiabatically 
along the contour C to give the real wavefunction \^f' P ) 
at point P' . Clearly this is an cigenstate of the Bloch 
Hamiltonian at this point in the Brillouin zone. A sep- 
arate way to obtain the eigenstate at P' is to apply the 
inversion operation on the state at P; I\^p). Again, 
the inversion operation can be constructed to yield a real 
wavefunction. Hence, these two wavefunctions can at 
most differ by a sign, 

1*$,) = o-(C)I\* P ) (Al) 

which is the inversion parity <r(C) of the curve C . Note 
that since I 2 = 1 <r(C) is independent of the direction 
of C. Although it depends on the band index, the band 
label is suppressed for clarity. 

Now consider two nodes N and N' being brought to- 
gether at M and choose two points P and P' along 
the perpendicular direction of this path. As shown in 
Fig|7]we choose two paths C\ and C2 connecting P and 
P'. Initially, C\ + C2 encloses a single Dirac node at 
N, and as a result we must have a(Ci)a(C2) = —1 be- 
cause the wavefunction must wind by n around a Dirac 
point. Now assume that the nodes annihilate on being 
brought together. Now, C\ + C2 encloses no singular- 
ity, so, the inversion parities c'(C) after this operation 
satisfy cr' {C\)a' {C2) = +1- However, since the wave- 
functions along C\ evolve smoothly during the annihila- 
tion we have: c'(Ci) = a(C\). Therefore we must have 




FIG. 7: We show the paths C\ and C2 connecting two points 
P and P which is inversion image of each other: (a) before 
annihilation and (b) after annihilation. C2 crosses the in- 
version symmetric point M, and before annihilation C\ + C2 
encloses one single Dirac node N. 

o~'(C2) = —<j{C2)- Note, by shrinking the curve C2 we 
can approach the M. point. Then, the inversion parity of 
the curve simply becomes the eigenvalue under inversion 
of the wavefunction at the Ai point: (Ji, 12). Therefore 
we conclude that in the node annihilation process, the 
inversion eigenvalue of each of the two states at the M. 
point changes sign. This is only possible if they have op- 
posite signs to begin with, I\l2 = — 1. In that case they 
can simply pass through each other, and the net result 
will be a sign change of the inversion eigenvalue of the 
higher and lower energy states. However, if they both 
have the same sign, = +1, it is not possible to af- 
fect a sign change. In this case, our assumption that the 
nodes annihilate is invalid - in fact a pair of bands with 
quadratic dispersion will touch at the M. point. 

If the two bands have opposite inversion eigenvalues, 
then the inversion matrix in the two bands is r 3 . And the 
inversion symmetric real Hamiltonian around M must be 
able to expand as (aSkx)T 1 + (e + bSk 2 , + c5k 2 , + d6k x Sk y )T 3 
to the quadratic order after choosing the k x axis to be 
along direction connecting the two nodes. We immedi- 
ately see that depending on the sign of e, the Hamilto- 
nian either has two or zero band touching nodes. This 
indicates that if the two bands have opposite inversion 
eigenvalues, the two nodes can always annihilate at A4. 

Finally we note that the SDW state in the five band 
model meets the conditions required for the above analy- 
sis to hold. Inversion T is a symmetry of the system, and 
the role of time reversal is played by T1Z' = SlZ(n — > 
—ft) o T7?.(defined in text) symmetric system, where n 
is the direction of collinear SDW. Let us choose the or- 
bital basis d a (a = 1 • ■ ■ n) to be eigenfunctions of inver- 
sion and label the eigenvalues to be I a . We then define 
the I o T1Z even basis d a in the following fashion: if 
I a = 1 then d a = d a , and if I a = —1 then d a = id a - 
For a collinear SDW with T1Z' symmetry the Hamilto- 
nian in d a is purely real, and so are the eigenfunction 
in the momentum space. In this basis, the arguments 
presented above can be made, leading to the conclusion 
that a pair of Dirac nodes cannot be annihilated at the T 
point. Hence, since we begin with six nodes in all, there 
will always be a leftover pair that is stable. 
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